Model-independent constraints on r 3 extra-interactions from orbital motions 
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Constraints on long-range power-law modifications [/pert oc r of the usual Newtonian gravita- 
tional potential Un oc r _1 are inferred from orbital motions of well known artificial and natural 
bodies. They can be interpreted in terms of a characteristic length £ which may be identified with, 
e.g., the anti-de Sitter (AdS) radius of curvature £ in the Randall-Sundrum (RS) braneworld model, 
although this not a mandatory choice. Our bounds, complementary to those from tabletop labora- 
tory experiments, do not rely upon more or less speculative and untested theoretical assumptions, 
contrary to other long-range RS tests proposed in astrophysical scenarios in which many of the phe- 
nomena adopted may depend on the system's composition, formation and dynamical history as well. 
Independently of the interpretation of £, the perihelion precession of Mercury and its radiotechnical 
ranging from the Earth yield £ < 10 — 50 km. Tighter bounds come from the perigee precession of 
the Moon, from which it can be inferred £ < 500 — 700 m. The best constraints (£ < 5 m) come from 
the Satellite-to-Satellite Tracking (SST) range of the GRACE A/B spacecrafts orbiting the Earth: 
proposed follow-on of such a mission, implying a sub-nm s -1 range-rate accuracy, may constrain £ 
at ~ 10 cm level. Weaker constraints come from the double pulsar system (£ < 80 — 100 km) and 
from the main sequence star S2 orbiting the compact object in Sgr A* {£ < 6.2 — 8.8 AU). Such 
bounds on the length £, which must not necessarily be identified with the AdS radius of curvature 
of the RS model, naturally translate into constraints on an, e.g., universal coupling parameter K, of 
the r" 3 interaction. GRACE yields K, < 1 x 10 16 m 5 s" 2 . 

PACS numbers: 04.50.-h, 04.80.Cc, 95.10.Ce, 95.10.Km, 96.30.-t, 91.10.Sp 



I. INTRODUCTION 



In this paper, we deal with power-law modifications of 
the usual inverse-square law [l|, Q . 

Extra dimensions, arising in string theory, supcrgrav- 
ity, M-theory Q and string inspired higher dimensional 
theories such as the brane- world models have re- 

cently gained increasing importance in physics in the 
context of the search for a quantum-gravity theory. 
In particular, large, non-compactified extra dimensions 
could potentially solve the long-lasting hierarchy problem 
0, S 0) II] • Indeed, if the standard model of particles and 
fields is restricted only on a (3+l)-dimensional brane, 
whereas gravity is allowed to propagate in the higher- 
dimensional bulk, the effective Planck scale in the four- 
dimensional spacetime can be made significantly larger 
than the electroweak scale, matching the experimental 
requirements. 

In the braneworld model by Randall and Sundrum (RS 
hereafter) our usual four-dimensional spacetime is a 
brane, which the standard model fields are constrained 
to, embedded in a five-dimensional anti-de Sitter (AdS) 
spacetime. In it, the fifth spatial dimension can be in- 
finite, with an AdS curvature scale I. Indeed, the RS 
model circumvents the need of compactifying all but the 
three observed spatial dimensions by including a bound 



state of the massless graviton on the brane @ resulting 
from the curvature, rather than the size, of the extra di- 
mension. At distances r S> i, the RS model implies a 
correction £/rs to the Newtonian gravitational potential 
of a body of mass M at second post-Newtonian order 
(2PN). It is given by 



U RS = —k 



GM ft 



(1) 
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where G is the Newtonian constant of gravitation, and k 
can assume different values depending on the schemes of 
regularization adopte d M . E.g., it can be k = 1 [f|, k = 
1/2 [Io[ , and k = 2/3 [III ]. The occurrence of a correction 
to the Newtonian potential of the form of Eq. ([TJ, which, 
however, is not necessarily limited to the RS model, being 
common to a wide class of power-law interactions 0], is 
important since, although it is only gravity that feels 
the presence of the extra dimensions, Eq. (jTJ allows for 
detectable effects on our (3+l)-brane that can be used in 
constraining the properties of the bulk. For upper limits 
on the brane parameter of other braneworld models, see, 

e.g., 0, m. 

Several large-scale tests of the RS model Q have been 
proposed so far in astrophysical scenarios; they claimed 
bounds on I at ~ 1 — 10 /im level, which is the same or- 
der of magnitude reached in laboratory-scale experiments 
fl4l [l5j . Anyway, such constraints typically depend on 
the particular interpretation of astrophysical observa- 
tions, and suffer from large systematic effects whose accu- 
rate knowledge is often lacking. Moreover, they arc often 
quite model-dependent in the sense that they heavily rely 
upon theoretical assumptions which are still speculative 
since they have not yet been tested independently with 
a variety of different phenomena, or have not yet been 
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directly tested at all. In our opinion, it is true also for 
those tests EMI 

requiring the least amount of informa- 
tion like, e.g., the evaporation of black holes Indeed, 
they are based on the application of a concept like the an- 
tide Sitter space/conformal field theory (AdS/CFT) du- 
ality (2ll | in braneworld gravity models that are asymp- 
totically AdS (such as the RS models). Moreover, the 
consequent theoretical prediction for the black hole evap- 
oration time (22T - [2~i| depends only on assumptions re- 
garding the braneworld model and, in particular, on the 
validity and implementation of the AdS / CFT correspon- 
dence in view of the "no-hair" conjecture [lHHfjjlj which 
is itself speculative and still awaits independent observa- 
tional checks [13, HH . Suffice it to say that the basis of 
the previously cited calculation for the black hole evapo- 
ration time [22h24| has been recently challenged in Ref. 
p9j . In addition to such theoretical considerations at 
fundamental level, it must also be remarked that the as- 
trophysical phenomena themselves used to constrain £, 
like the orbital evolution of black-hole X-Ray binaries or 
the behavior of black holes in extragalctic clusters, are 
not lacking of uncertainties, both from a theoretical and 
observational point of view. E.g., they may crucially de- 
pend on the composition, formation and dynamical his- 
tory of the systems considered. Last but not least, black 
holes may well not exist at all [3(J EJ; e.g., even in the 
case of the compact object in Sgr A*, there are not yet 
direct evidences that its ~ 10 6 M Q mass is actually con- 
centrated within its Schwarzschild radius R s = 0.084 AU 
[3^ . A signature for the absence of event horizons was 
even looked for by the authors of Ref. [33| . In conclusion, 
to date, a definite proof for the existence of Kerr black 
holes is still lacking despite a wealth of observational ev- 
idence 34j. Moreover, from a broader point of view, in 



correctly assessing the relevance of the laboratory-based 
tests of long-range modified models of g ravity it should 
be mentioned that some authors [35l [36| pointed out the 
need of go to space as well since such extra gravitational 
degrees of freedom may have an environmental depen- 
dence. 

Thus, in Sec. [TT] we will use well known and largely 
tested orbital motions (see Sec. Ill A|) of some natural 
(Sec. Ill Bp and artificial (Sec. Ill C|) bodies in the solar 
system to infer constraints on the AdS radius of curvature 
which, if on the one hand, are not at /jm level, on the 
other hand can certainly be considered less speculative 
than those obtained in astrophysical contexts. Indeed, 
apart from the fact that there are no doubts about the ex- 
istence and the properties of the planets of the solar sys- 
tem and of man-made Earth's satellites, the dynamical 
effects used to constrain I are straightforwardly obtained 
from Eq. ([1]) , without any additional hypothesis concern- 
ing untested phenomena. Moreover, competing effects 
acting as systematic errors are known with a compara- 
tively much better accuracy. We notice that McWilliams 
[37| recently proposed to constrain £ via gravitational 
wave measurements in the solar system: LISA would al- 
low to place bounds on the AdS radius of curvature of 



the order of £ ~ I /im from the event rate of stellar black 
holes inspiraling gravitationally into supermassive black 
holes, and of £ < 5 /zm from the observation of individ- 
ual galactic binaries containing a stellar mass black hole. 
In Sec. Ill Dl we will also use the well known, and ex- 
tensively studied, double pulsar binary system and the 
main sequence S2 star orbiting the compact object in 
Sgr A*. Finally, we stress that our results are not nec- 
essarily limited to the RS model @, being valid also for 
other theoretical schemes yielding power-law interactions 

cx r- 3 [IllalSi. 



II. RS LONG-TERM ORBITAL EFFECTS AND 
COMPARISON WITH THE OBSERVATIONS 

A. Analytical calculation of the secular precession 
of the pericenter 

The long-period effects caused by Eq. ([T]) on the or- 
bital motion of a test particle of mass m can be computed 
perturbatively by adopting the Lagrange equations for 
the variation of the osculating Keplerian elements [43| : 
their validity has been confirmed in a variety of indepen- 
dent phenomena. Generally speaking, they imply the use 
of a perturbing function 1Z which is the correction £/p er t 
to the standard Newtonian monopolc term. In the case 
Expert = Urs, the average over one orbital revolution of 
the perturbing function 1Z is straightforwardly obtained 
by using the true anomaly / as fast variable of integra- 
tion: it is 



(n) = - 



GMkt 2 



a 3 (I - e 2 ) 



2^/2' 



(2) 



where a is the semimajor axis and e is the eccentricity of 
the test particle's orbit. From Eq. ([2]) and the Lagrange 
equation for variation of the longitude of pericenter vo 



dm 
~dt 



1 



1 



,V2 



d{K) tan (J/2) d{K) \ 
de + (l- e 2)V2 dl J ' 

(3) 



it turns out that w experiences a secular precession given 
by 



dvo 
~dt 



3(GM) 1/2 k£ 
a 7/2 (! _ e 2y 



(4) 



In Eq. ©, n b = ^GAI/a 3 is the Keplerian mean motion 
and I is the inclination of the orbital plane to the refer- 
ence {2, y} plane. The longitude of pericenter vj = Q + lu 
is a "dogleg" angle since it is the sum of the longitude of 
the ascending node f2, which is an angle in the reference 
{x, y} plane from a reference x direction to the intersec- 
tion of the orbital plane with the {x, y} plane itself (the 
line of the nodes), and of the argument of pericenter w, 
which is an angle counted in the orbital plane from the 
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line of the nodes to the point of closest approach, usually 
dubbed pericenter. The precession of Eq. Qj, which is 
an exact result in the sense that no a-priori assumptions 
on e were made, agrees with the one obtained by Adkins 
and McDonnel in Ref. [45| with a more cumbersome cal- 
culation. To facilitate a comparison between such two 



results, we note that, in general, the authors of Ref. [45 1 
work out the perihelion advance per orbit A6 p : it corre- 
sponds to (w) Ph-, where Pb = 27r/n^ is the orbital period. 
Moreover, in the potential energy V(r) = a_(.,- +1 )r~^ +1 ) 
it must be posed j = 2 and a_3 — > —GMmki 2 , while 
in A0 p (-(j + 1)) of Eq. (38) in Ref. (H it must be 



set L —> a(l — e ), X2( e ) = 6. With such replacements, 
it can be shown that the advance per orbit of Eq. (38) 
in Ref. [45j corresponds just to our precession in Eq. 
(U]). Moreover, it turns out that the analytical result of 
Eq. ([4| is confirmed by a numerical integration of the 
equations of motion for Mercury with, say, £ = 10~ 6 AU 
and fc = 1/2: both yield 5.4 milliarcseconds per century 
(mas cty - 1 hereafter) . The choice of the numerical value 
adopted for £ is purely arbitrary, being motivated only 
by the need of dealing with relatively small numbers. 



B. Constraints from solar system planetary orbital 
motions 

The corrections Avu to the standard Ncwtonian- 
Einsteinian secular precessions of the longitudes of the 
perihelia are routinely used by independent teams of as- 
tronomers [46l |47| as a quantitative measure of the maxi- 
mum size of any putative anomalous effect allowed by the 
currently adopted mathematical models of the standard 
solar system dynamics fitted to the available planetary 
observations. Thus, Art can be used to put constraints 
on the parameters like I entering the exotic models one is 
interested in. From Eq. (01 it turns out that the tightest 
constraints come from Mercury, which is the innermost 
planet with a — 0.38 AU, for k = 1. 

Ficnga et al. [47[, who used also a few data from the 



three flybys of MESSENGER in 2008-2009, released an 
uncertainty of 0.6 mas cty -1 for the perihelion precession 
of Mercury, so that it is £ < 34 - 48 km for k = 1 - 1/2. 
The uncertainty in the pre-MESSENGER Mercury's per- 
ihelion extra-precession by Pitjeva [46| is about one order 
of magnitude larger (5 mas cty -1 ). 

Slightly tighter constraints on £ come from the in- 
terplanetary Earth-Mercury ranging. Indeed, according 
to Table 1 of Ref. [47], the standard deviation a^ p of 
the Mercury range residuals Ap, obtained with the IN- 
POPlOa ephemcrides and including also three Mercury 
MESSENGER flybys in 2008-2009, is as large as 1.9 m. 
A numerical integration of the equations of motion of 
Mercury and the Earth yields a RS range signal with the 
same standard deviation for £ < 13 — 18 km (fc = 1 — 1 /2) . 
Figure [1] displays the case fc = 1/2. It turns out that the 
residuals of right ascension (RA) and declination (DEC) 
in Ref. [47| yield much less tight constraints. 



As far as natural bodies of the solar system are con- 
cerned, the Moon yields better results. Its orbit is ac- 
curately reconstructed with the Lunar Laser Ranging 
(LLR) technique [H since 1969; Figure B-l of Ref. 
[491 ] shows that the residuals of the Earth-Moon range 
are at a cm-level since about 1990. The secular pre- 
cession of the lunar perigee is known with an accuracy 
of about 0.1 mas yr _1 [50l453j . so that Eq. (@| yields 
£ < 524- 741 m (ft = 1- 1/2). 



C. Constraints from the GRACE spacecraft 
orbiting the Earth 

Remaining within the solar system, tighter constraints 
can be obtained from selected spacecrafts orbiting the 
Earth. The Gravit y R ecovery and Climate Experiment 
(GRACE) mission [HJ], jointly launched in March 2002 
by NASA and the German Space Agency (DLR) to map 
the terrestrial gravitational field with an unprecedented 
accuracy, consists of a tandem of two spacecrafts mov- 
ing along low-altitude, nearly polar orbits continuously 
linked by a Satellite to Satellite Tracking (SST) mi- 
crowave K-band ranging (KBR) system accurate to bet- 
ter than 10 um (biased range p) [551 ] and 1 p,m s _1 (range- 
rate p) [55|,l56[. Studies for a follow-on of GRACE show 
that the use of a intcrfermotcric laser ranging system may 
push the accuracy in the range-rate to a ~ 0.6 nm s -1 
level. A numerical integration of the equations of motion 
for GRACE A/B, including also the mismodelled signal 
of the first nine zonal harmonics of geopotential [43[ ac- 
cording to the global Earth gravity model GOCO01S [58j |. 
shows that the SST range is more effective than the SST 
range-rate in constraining the RS parameter for which 
it holds I < 5 m. Figure [2] depicts the numerically in- 
tegrated GRACE SST range signal due to Eq. Q and 
the aforementioned mismodeled zonals. It can be shown 
that a sub-nm s _1 level of accuracy in the GRACE SST 
range-rate would imply the possibility of constraining £ 
down to ~ 10 cm level. 



D. Constraints from the double pulsar and Sgr A* 

Constraints comparable with the planetary ones can be 
obtained from the periastron of the double pulsar PSR 
J0737-3039A/B system [H, Hi Indeed, the semima- 
jor axis of its relative orbit amounts to just a = 0.006 
AU (olj . The present-day accuracy in measuring the 
secular precession of the periastron is 6.8 x 10 -4 de- 
gree per year (deg yr" 1 in the following) 6l[; thus, a 
straightforward application of it to Eq. (0| would give 
£ < 11 - 16 km (fc = 1 - 1/2). Actually, the larger 
uncertainty in the theoretical expression of the general 
relativistic 1PN periastron precession must be taken into 
account as well. It is as large as 0.03 deg yr _1 [g2], so 
that it yields £ < 80 - 112 km (fc = 1 - 1/2). 

The pcrinigricon of the S2 star, orbiting in 15.98 yr the 
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Earth-Mercury RS range perturbation (L = 18 km, k = 1/2) 




t(yr) 



FIG. 1. 2008-2009 Earth- Mercury range shift Ap, in m, due to an unmodeled potential like Eq. fTJ) for t = 18 km and k = 1/2 
computed from the difference between two numerical integrations of the equations of motion in cartesian coordinates performed 
with and with out Eq. ([!]). Both the integrations s hare the same initial conditions retrieved from the NASA WEB interface 
HORIZONS at|http://ssd.jpl.nasa.gov/horizons.cgi| for the epoch J2000.0. It is (Ap) = -0.04 m, a Ap = 1.9 m. 



GRACE RS SST range perturbation (L = 5 m, k = 1/2) 




-1 

-2 
1-3 
*-4 
-5 
-6 




V 



0.2 0.4 0.6 0.8 1 

t(d) 



FIG. 2. Daily GRACE SST range shift Ap, in pm, due to an unmodeled potential like Eq. JTJ (£ = 5 m. k = 1/2) 
and the first nine mismodelled zonals of geopotential according to GOCO01S [SH] computed from the difference between 
two numerical integrations of the equations of motion in cartesian coordinates performed with and without the pertur- 
bations considered. Both the integrations share the same initial conditions contained in the files GNVlB_2003-09-14_A_00 
and GNVlB_2003-09-14_B_00 retrieved from |ftp: / / cddis.gsfc.nasa.gov/pub/slr/predicts/current / graceA_irvs_0 81202_0.gfz| and 
ftp: / /cddis. gsfc.nasa.gov / pub /sir / predicts/current /graceB_irvs_081201_l.gfz. The epoch is 13 September 2003. (See 
ftp://podaac.jpl.nasa.gov/pub/grace/doc/Handbook_lB_vL3.pdf for the explanation of the GPS Navigation Data Format 
Record (GNV1B) format). It is (Ap) = —3 pm, ctap = 2 pm. 



compact object hosted in Sgr A* with M = 4x 10 6 M Q 
[fiU at (r) = 1433 AU from it J6_|, yields much weaker 
constraints on the asymptotic AdS curvature. Indeed, 
an accuracy of [64[ 0.05 deg yr -1 can be inferred for the 
pcrinigricon precession of S2: thus, I < 6.2 — 8.8 AU. 



III. MODEL-INDEPENDENT CONSTRAINTS 
ON r" 3 INTERACTIONS 

The validity of the previous results is not necessarily 



limited just to the RS braneworld model. Indeed, they 
are, in fact, quite mo del- independent in the sense that 
they hold for any long-range modification C/ por t oc r~ 3 
[3, [38-42] of the usual Newtonian potential. In partic- 
ular, t must not necessarily be identified with the AdS 
radius of curvature of the RS braneworld model. 



By assuming that the putative, new interaction does 
not depend on the specific matter distribution generat- 
ing the gravitational field, a universal parameter can be 
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introduced with the replacement 

kGMf -> K, [K\ = L 5 T~ 2 . (5) 

Thus, the bounds on £ of Section III EH Scction III Dl 
straightforwardly translate into constraints on tC itself. 

It turns out that the tightest bounds come from the 
Earth-GRACE system yielding K < 1 x 10 16 m 5 s" 2 . 



IV. SUMMARY AND CONCLUSIONS 

Table U and Table [TTI resume our findings. 

As far as the £ 2 form of the r~ 3 extra-potential is con- 
cerned, Table Q] tells us that tightest bounds, at the m 
level, come from the Satellite-to-Satellite Tracking rang- 
ing between the GRACE spacecrafts. 

If we interpret them in terms of the RS model, from 
a general point of view, we find not entirely adequate 
arguing that Earth-based laboratory-scale tests of sev- 
eral modified models of gravity are in all respects su- 
perior since they can deliver much tighter constraints. 
Indeed, apart from the basic fact that it is important 
to scrutiny a theoretical paradigm in different, comple- 
mentary scenarios, some authors also pointed out that 
there are theoretical reasons for wanting to test such 



foundational issues of gravity in space since extra gravi- 
tational degrees of freedom could have an environmental 
dependence. Moreover, about the seemingly superior- 
ity of certain constraints of the RS model inferred from 
some astrophysical systems, often they rely upon more 
or less speculative and untested theoretical assumptions, 
and the phenomena adopted may depend on the system's 
composition, formation and dynamical history as well. 
Thus, we reiterate the importance of finding solid, well- 
understood means of constraining £ like those proposed 
here. We also point out that the same method introduced 
here could yield another factor of 10 2 improvement, con- 
straining £ < 10 cm, provided that the proposed improve- 
ment to satellite-to-satellite tracking is implemented for 
the follow-on of the GRACE mission. 

Last but not least, we remark that, actually, our anal- 
ysis is not necessarily limited to the RS model since it 
is valid for whatsoever theoretical scheme predicting r~ 3 
corrections to the Newtonian potential. In this respect, 
the occurrence of a characteristic length scale £ is com- 
mon to a wide class of power-law interactions, so that 
it should not necessarily be thought of as the AdS ra- 
dius of curvature in the RS model. Thus, the bounds 
of Table U have a wider range of applicability. Table 
HT1 shows that GRACE yields the tightest bounds on a 
universal coupling parameter K, which is constrained to 
< 10 16 m 5 s~ 2 . It could be pushed down to a 10 12 m 5 s~ 2 
level by GRACE follow-on. 
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